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Abstract 

The effect of the contribution of the vertical part of the real time path is studied completely 
, in the case of two points functions. Indeed, this vertical part generally contributes in the cal- 

Qj I culation of a given graph. Moreover, this contribution is essential in order to have a consistent 

i-S^ ' equilibrium theory: thanks to this contribution, the Green function are effectively invariant 

by time translation, as they should be. As a by product, it is shown that the perturbative 
^ I calculations give a result which does not depend on the initial time ti and final time tp of the 

■ path. The property of independence with respect to tj is closely related to the KMS condi- 

^ tions, i.e. to the fact the system is in thermal equilibrium. In the case of two point functions, 

the contribution of the vertical part can be taken into account by the n(|/co|) prescription in 
the usual RTF Feynman rules. 
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I. Introduction 



Quantum Field Theory can be formulated on every descending path in complex time plane, 
provided it begins at an initial time ti and ends at the time ti — i/3. Among these paths, of 
particularly important use are the so-called real-time paths, represented on figure 1. This kind 
of path has two parts parallel to the real time axis, and two vertical parts. In the following, 
we will only consider in this familly of paths the path corresponding to a = 0*^^ (see figure 2), 
so that there will be only one vertical segment, which will be refered in the following as the 
"vertical part" and denoted simply Cy instead of C4. 

It is necessary to find a proper way to get rid of the vertical part if one wants to have simple 
Feynman rules in Fourier space. The usual way to achieve that [1 ,2 ,3 ,4 ,5 ,6 ,7] has been to 
simply drop the vertical part, arguing that it does not contribute to the Green functions in the 
limit tj —00, tp +00, thanks to the Riemann-Lebesgue lemma which is said to imply the 
factorization of the generating functional: Z[j] = Zi2[j] Zzi[j]. However, it has been noticed 
that the vertical part can contribute to the Green functions, for example in vacuum graphs 
[8 ,9] and even in every two point function [10]. Alternative paths with only two branchs have 
been introduced [11] in order to justify the 2x2 matrix formalism, but the limiting procedure 
to obtain Real Time Formalism in Fourier space (simply called RTF in the following) is not 
completely transparent. In section II, we see by a canonical approach (i.e. without the path 
integral formalism) how this vertical part does arise, and what is the precise meaning of its 
suppression. We show that dropping the vertical part is equivalent to perform averages with 
the free density operator exp(— /3iJo) instead of exp{—(3H), the dynamics of the fields still 
being governed by the complete Hamiltonian. 

Another problem related to the RTF is the so called n(|A;o|) prescription. Indeed, the free 
propagators contain the statistical function n(-) always accompanied with a factor S{kQ — 
(jj^), where ojk = Vk'^ + m"^- So, in the free propagators, it makes no difi:erence to have 
n(|fco|) or n{ujk)- But, because the cancellation of pathologies generates derivatives of the 
5{-) distributions, the choice of the argument of the statistical factors is very important in 
perturbative calculations where pathologies arise. It has been shown that in order to preserve 
the KMS relations when the propagators are regularized by a finite parameter e > 0, the 
n(]/col) prescription is necessary [3,11]. 

In [10], it is noticed that the vertical part contributes to every connected two point diagram, 
and that this contribution can be included among the other terms (those which come from the 
horizontal parts of the path) by simply substituing in them n(|/i;o|) to n[uJk)- So it seems that 
the two announced problems are closely related, and that the n{\ko\) prescription is in fact 
the consequence of the vertical part contribution. However, the proof given in [10] in order 
to establish this property for the general two point diagram still uses the limit ti — > —00, 
tp — ^ +00 in association with the Riemann-Lebesgue lemma. 

In section III of the present paper, we perform simple calculations on the complete real 
time path, with finite ti and tp, and then Fourier transform the result in order to see what is 
the effect of the vertical part. We reobtain the result announced by [10]. Moreover, we also 
see that the result of the calculation is independent of the times tj and tp, so that the proofs 
based on the limits tj —00, tp +00 seems suspect. In section IV, we prove explicitly that 



Nevertheless, most of the calculations made in this paper are valid whatever is the path 
going from ti to tj — i(3. 
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every diagram is independent of tj and tp- The independence with respect to tp is simply 
related to causality, because tp must be chosen greater than any other time in the diagram. 
The independence with respect to tj is related to the fact that the system is in equilibrium, 
and is proven by use of the KMS relations. 

In section V, we prove that the vertical part contributions can be taken into account by 
the n(|A;o|) prescription for every two point function, without invoking the ti — > — oo and 
tp ^ +00 limits. 

In section VI, we show that it is possible to have a contribution of the vertical part in n 
point functions (n > 2), for certain configurations of the external momenta. 



II. Canonical derivation of the generating functional 



The path-integral technique, in association with the Feynman-Matthews-Salam formula [3], 
is the most powerful way to quantize a classical field. Nevertheless, in the case of thermal field 
theories, a canonical approach can clarify the actual role of the vertical branch of the path. 
We reproduce here this derivation in the simple case of a real scalar field^^' [3 ,12 ,13] . 

2.1 The generating functional for a general average 

Let (f){x) be the field operator in the Heisenberg picture, and (l)i{x) the field operator in the 
interaction picture, the interaction being switched on at the time ti. They are related by: 

(cl)ix) = U{ti,t)cl>,ix)Uit,ti) 



t2 



X 



[1) 



where the integration is performed along a path going from ti to t2 in the complex-time plane, 
and where £/(•) is the interacting part of the Lagrangian. The operator U{ti,t2) verifies the 
following basic properties: 

\U{h,t2)U{t2,h) = U{tl,h) 



By introducing an arbitrary time tp, we can write: 

0(x) = U{ti,tp)U{tp,t)(j),{x)U{t,ti) 

from which we deduce: 



Tc(t){xi)---(j){Xn) = Tc 



4>j{xi) ■ ■ ■ (l)j{xn) expi J d^x Ci{(l)i{x)) 



(3) 



(4) 



Of course, the case of more complicated fields like non abelian gauge fields is prohibitively 
tedious to treat by this method. But we expect that the interpretation of the vertical part 
obtained here for scalar fields still holds for other kind of fields. 
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where Ci is a path going along the real time axis from tj to tp, and C2 the reverse path from 
tp to tj. 

Now, we define the Green functions to be the averages of the path ordered products calcu- 
lated in (4), an average operator being in general a linear form a;(-), normalized to a;(i) = 1. 
Prom (4), we can directly write the generating functional of these Green functions: 



1 _J_ 

Z^[0] iSj{xi) iSj{xr 



u;(T,0(xi) • • • ct>ixn)) = " " " -TTTZr^Z^ i^h^o (5) 



Z'^ij] = u Tecxpi J d^x CiiMx))+jix)(p,ix) (6) 
In particular, we will denote: 

where Hq is the free part of the Hamiltonian H = Hq + H' in the Heisenberg picture. 

2.2 The introduction of the vertical part 

The fundamental formula [13], useful to perform equilibrium averages (i.e. with Uh{-)), is: 

^-PH ^ e-^^°Tcexpi J (fx Ci{(t),{x)) (8) 

where is a path going from tj to tj — ifH, and which will be taken parallel to the imaginary 
time axis (figure 2). Then, thanks to the formula (8), we can rewrite the generating functional 
of the equilibrium Green functions as: 

Z^^ij] = ^^J"h] ^h, ^Tecxpi j d^x Ci{cP,{x))+j{x)(t>,{x)^ (9) 

where C = Ci © C2 © is the real time path with cr = 0. If, instead of calculating equilibrium 
averages, we decide to perform ujuoi ) averages, we do not need to use the formula (8) and we 
can write directly: 



=a;^ jTeexpi j 



d^xCM,{x))+j{x)<t),{x)\ (10) 



So, apart from a ^-independent factor which plays no role in the perturbative calculation 
of a Green function with (5), we see that the only difi:erence between the last two formulae 
is the integration path. The last one shows that when we neglect the vertical part, we are in 
reality calculating averages with the free density operator exp{—(3HQ). 
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We can further simplify these generating functionals by the same manipulations as the 
zero temperature ones [14], and by using the Wick's theorem which remains valid for iOHf^{-) 
averages of free fields (see appendix 1). Indeed, we have: 

Uh^ (Xcexpi J d^x j{x)(j)i{x)\ = exp-^J d^x d^y j{x)j{y)uH^{Tc(j)j{x)(j)j{y)) (11) 

\ 7 / 7 

where 7 is an arbitrary path, so (11) applies both to (9) and (10). (11) implies that to evaluate 
a graph, we need to attach a free propagator Go{x,y) = u!Ho{T^c4>i{x)(f)i{y)) to each line, and 
to integrate over the coordinates of each vertex, the integral over time being performed along 

the path 7. 

So, the important conclusion of this part is that the suppression of the vertical part of 
the time path is equivalent to perform only ujho{-) averages. More precisely, we can say that 
adding the vertical part allows us to do at the same time the perturbative expansion in the 
dynamics of the fields (in the same way as in the zero temperature expansion) and in their 
statistics (i.e. in the density operator, which also contain the coupling constant): the vertical 
part enables us to have a true equilibrium theory at each order of the perturbative expansion. 



III. A simple example: the insertion of mass terms 



3.1 Calculation in time coordinates 
3.1.1 Introduction and notations 

In order to see what the effect of the vertical part in effective calculations is, we will do the 
simplest calculation which consists in the convolution of two propagators: 

F{xi,X2)= J d'^x Go{xi,x)Go{x,X2) (12) 

c 

This is typically the calculation one would do to insert a mass term on the propagator of a 
particle. We can further simplify the computations by using a Fourier representation with 
respect to the spatial dependence: 



G 

It is well known that'^': 



o{ti,t2;k) = J d^xe-"'-'' Go{ti,0,t2,x) (13) 



Go{h,t2;ujk) ^ Go{h,t2;k) = ^\e-^'^''^*^-*'\ec{t2 - h) + nB{ujk))+ 

e-'^'^''-''\0c{h-t2)+nB{u;k)) 



Whenever a two point function depends on k only through ujk, we will recall this by 
replacing k by tUk in the list of its arguments. 
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So, the Feynman rules will be slightly modified and we have now to compute: 



F{ti,t2;u;k) = J dtGo{ti,t;ujk)Go{t,t2;ujk) (15) 
c 

3.1.2 The tF-independence 

Before going to the calculation itself, we make a general remark on the independence of 
the Green functions with respect to the time tp- To see this property, we can choose another 
final time t'p > tp- Then if Gtp is the Green function calculated on a path which final time is 
ip, we will have: 



Gt'^ - Gtp 



= J g{t)dt (16) 



where 7 is a path going from tp to t'p and then back to tp and g{t) a function which does not 
contain explicitly tp or t'p. Moreover, if tp is greater than all the other times in the problem, 
g{t) will take the same values on the upper and lower part of 7, and will be equal to Gtp, 
as announced. 

3.1.3 Result and commentaries 

After some straightforward calculations, we find: 

F{t,,t2;iVk) = { nlM [ ^(e-'^^*-*^) + e--'^^*-*^)) 
+ riBiuJk) [ {t2 - ti)(e-^'^'=(*^-*^) - e^'^'=(*^-*^)) 

_|_ ^1 r^-iuJkit2-ti) _^ ^iuJk{t2-tl)\ 



_J_^^ioJk{ti+t2-2tj)^-2uJkl3 _j_ g-iwfe(ti+t2-2tj)^ 



+ [ 0cit2 - ti){t2 -ti + ^)e-^'^'=(*--*0 + 0^(^ti - t2){ti -t2 + — )e^'^'=(*--*^) 
2iujk 

Firstly, we notice that there appears some terms which are not invariant by time translation, 
depending on ti + t2 — 2ti. But, we verify that the sum of all these terms is zero. Therefore, 
the final result is invariant by time translation and independent of the initial time tj. We 
must precise that the contribution of the vertical part*^^' (terms in and in exp(it2/3u;fc)) is 



The separation between the contributions of the horizontal part and of the vertical part 
has no intrinsic meaning: it depends on the position of the external times on the path (if they 
are all on the vertical part, the entire result comes from the vertical part). So, here and in 
the following, the contribution of the vertical part is defined to be the one obtained when the 
external times are on the horizontal branchs. 
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essential in order to have this cancellation. We will generalize this property to a general graph 
in the next section. Another property of this result is the fact that it remains contributions 
of the vertical part (terms proportional to and we have no possibility to make them 

disappear by the limit procedure tj — > — oo, tp — > +00 because the result does not contain ti 
or tp anymore. 

We can also verify explicitly the mass-derivative formula [10 ,15] : 



i^^G'o(*i,*2;t^fe) = J dt Go{ti,t;ujk)Go{t,t2;u;k) (18) 

c 

It is important to notice here that this formula holds only if the integration is performed on 
the complete path C, because the left-hand side of (18) contains those terms in exhibited 
above. 



3.2 The result in Fourier space 

Now, the following step is to perform the Fourier transformation of F{t\,t2] k) with respect 
to the difference t2 — ti**^' in order to see which are the Feynman rules in energy space (they 
are not obvious because of the contribution of the vertical part in the result). We will denote 
Dij^ujik) the usual real-time matrix propagator (with a = 0) with distribution function nB{iOk) 
and -Dij,feo (^) the same propagators with distribution function nB{\ko\) (see appendix 2). These 
two kinds of free propagators are equal because each distribution function is accompanied 
by a 5{kQ — a;|) factor. Then, by a simple but tedious calculation, we find for the Fourier 
transformation : 



where the last term is the Fourier transformation of the vertical part contribution. Moreover, 
in the first term, the cancellation of pathologies makes derivatives of delta distributions appear 

(see appendix 3 and section V), so that this term will not be the same if we replace tOk by /cq 
in it. But, fortunately, if we do such a substitution, we obtain a supplementary term which 
cancels exactly the second term of (19), so that [10]: 



Fij{ko,Uk) = {DkMrzDk,{k)).. (20) 

In this simple example, we have seen that dropping the vertical part contribution is equi- 
valent to use type propagators in the usual real-time Feynman rules [1 ,2 ,3 ] , and that we 
can take this contribution into account by using fco~type propagators with the same Feynman 
rules. 



This Fourier transformation is legitime once we have verified that the result is invariant by 
time translation. Moreover, the function F{-) being independant of the times ti and tp, we 
are allowed to perform an usual Fourier transformation with times varying from — cx3 to +00. 
We then define a 2 x 2 matrix according to the position of the external times on Ci or C2, in 
the usual way. 
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3.3 Relation with part II 

To end this section, we verify the physical interpretation of section II, i.e. the fact that 

dropping the vertical part is equivalent to perform averages with the free statistics. For that 
purpose, we treat perturbatively an additional mass term —ji^cfp' /2 (figure 3). If we write the 
free matrix propagator as: 

n fi\ TT ( i{k'^ - m'^ + ie)-^ \„ 

where the matrix U contains the distribution functions and verifies Ut^U = t^. The only effect 
of the summation of all the extra mass insertions is to replace by + jj,^ in the central 
matrix. Consider for example the 11 component of the resummed propagator: 

'^iiMkolik) = fe2_^2'_^2 + ^, + '^^nsMkoimk' -m'- fi') (22) 

where ujk = Vfc^ + Tn"^ corresponds to the former mass. We then see that if we use the 
a;— type free propagators, the result corresponds to the new dynamics (through the change 
of pole) and to the former statistics (corresponding to the mass m), as announced. If, on 
the contrary, we use the /cq— type free propagators, the resultant statistical distribution will be 
nBi^/k"^ + m"^ + A*^), i-e. the result takes into account the change of mass also in the statistics. 

3.4 Conclusion 

On the simple example treated in this part, we have seen that the vertical part does 
contribute to the final result, and that taking into account this contribution is equivalent to 
the so-called n(|fco|) prescription. Moreover, we insist again on the fact that this contribution 
is independent of the times tj and tp- We will extend later this property to a general two 
point function. 



IV. A general proof of the tj-invariance 



4.1 Proof in the case of scalar fields 

In the preceding section, we noticed on an example that the result of the time integration 
along the complete path C was independent of ti, and was invariant by time translation. Before 

going on to a general proof of this result, we must emphasize the fact that this is in complete 
agreement with our hypothesis of thermal equilibrium. Indeed, a system in equilibrium must 
be invariant by time translation and can't have memory of its initial time: therefore its Green 
functions must verify also these two properties. We must also insist on their non obvious 
character: indeed, the free propagators verify them, but the integration path does not''^'. 



'''' As we will see later, we need an extra property of the propagators expliciting the fact that 
they are equilibrium propagators: i.e. the KMS relations. 
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In the following, we will denote the value of a n point arbitrary graph by: 



Gtj{tl, ■ ■ ■ ,tn) ^ / C^^n+i -P'(il) • • • 5 in! in+1) • • • 7 ^n+v) (23) 

where v is the number of vertices in the graph, and F(-) the product of the free propagators 
corresponding to the lines of the graph, before we perform the vertex-time integrations***' . The 
tj index is there to recall that the integrations are performed on a path beginning at tj and 
ending at tj — i(3. 

Firstly, we can prove the equivalence between the following two properties: 

V(ii,- Gtj{h,---,tn) = Gt,{ti + a,---,tn + a) (24) 

V(tl, • • • , tn), Gti{tl, ■ ■ ■ ,tn) = Gtj-aitl, ■ ■ ■ ,tn) (25) 

because the function F{-) depends only on time-differences (see (14)). Therefore, the invariance 
by time translation is equivalent to the independence with respect to the initial time tj. We 
now prove this last property, first for scalar fields and then for fermions, eventually with a 
chemical potential. 

For that purpose, we need the following two lemmas: 



d 
dt 



- j f{t)dt = f{tr-m-f{ti) (26) 



where / is a function defined on the path Ctj, which does not depend explicitly on tj, and: 

(Go{ti - i(3,t2;u)k) = Go(ti,t2;u;k) ^^j^ 
\Go{ti,ti - i^;ujk) = Go{ti,tr,u;k) 

The last one is only a consequence of the KMS relations, that is to say of the thermal equili- 
brium property of the system. 

Then, the F(-) function verify also this kind of KMS property: 

Vj G [1, • • ■ ,v], F{- ■ ■ ,tn+j-l,ti - il3,tn+j+l, ■ ■•) = F(- • • ,tn+j-l,ti,tn+j+l, • • •) (28) 

It is therefore easy to verify by recursion that the v integrations will lead to a function Gtj{-) 
which does not depend explicitly on the initial time tj. 



(8) We have not written the integrations over the three momenta of independent loops, because 
the spatial dependence of the Green functions has nothing to do with the properties we are 
looking at here. 
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4.2 Extension to fermions and chemical potential 

The case of fermions with chemical potential is a bit more intricate because the KMS 
relations are now: 

(Soiti,tr,iOk) = -e-^^'"'So{ti,ti -i(3-iOk) 
\So{ti,t2;u;k) = -ef^'"'So{ti-i(3,t2;u;k) 

where /j, is the chemical potential and q the conserved charge carried through the propagator 
from its first point to its second point. Because the number of fermionic lines going to a vertex 
is even, the minus signs cancel mutually when we try to establish (28) with fermions. In the 
same way, the exp(±/3//g) factors cancel because of the nullity of the total conserved charge 
arriving at a vertex. So, (28) is still true in the general case of fermions with chemical potential 
(of course, the KMS relations with respect to the external times depend on the type of the 
external fields), and that was the only property of the F(-) function we need to prove the 
invariance of Gtj{ ) with respect to t/. 

4.3 Application to ITF 

This paragraph lies a bit outside the general line of this article: it shows how the preceding 
arguments can lead to define "ITF-like" Feynman rules, not only on the usual imaginary time 
path, but also on every other path going from tj to tj — i(3. The property of t/— invariance 
can be used to introduce the discrete Matsubara frequencies of the imaginary time formalism. 
Let's assume that we have performed the calculation of a graph in time coordinates: we are 
led to a function G{ti, ■ ■ ■ ,tn) which does not depend on tj and which verify KMS relations 
according to the type of fields attached to the external lines. The ITF corresponding function 
can be defined as [16]: 

„ n 

G(a;i,---,a;„)= J {dU e'''''*) G{h, ■ ■ ■ ,tn) (30) 



i=l 



Naturally, we require that this ITF function should be independant of the initial time tj. It 
is this requirement that will lead us to the discretization of the Matsubara frequencies. Let's 
denote: 

F(ii,---,^,)= (TTe-^^O G{h,---,tn) (31) 




According to the proof given in the preceding paragraphs, the ITF function will be independent 
of ti if the function F(-) verify scalar-like KMS relations: 

Vj G [l,---,n], F{- ■ ■ ,tj-i,ti,tj+i,- ■ ■) ^ F{- ■ ■ ,tj-i,ti - i/3,tj+i,- ■ ■) (32) 

This lead us to the usual frequencies: 
(i) if the line n°i carries a real scalar: 

ooi = with Hi G Z (33) 
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(ii) if the line n°i carries a fermion: 



(2ni + l)7r . , „ 

cui = - — ^ with ni G Z (34) 

(iii) if the hne n°i carries a fermion with conserved charge qf. 

(2ni + l)7r . . , „ 

cjj = i/iQi With Hi E Z (35) 

It is then easy to prove the following property [16]: 

G(u;i, ■■■,ujn) = -i/35^,oG{uj2, • • • ,t^n) (36) 

with 

/n 
Ylidtie-'^*'*) G{0,h,---,tn) (37) 



i=2 



where S = cji + • • • + aj„ is the total energy entering the graph. We can note that the 
invariance of the result of time integrations with respect to tj is essential in this proof. This is 
quite normal because the t/— invariance is related to the invariance by time translation, which 
is itself related to the conservation of total energy we want to prove here. 
The inverse transformation is: 

G{h,---,tn)=(^) J2 [fle'^^A Giuji,---,un) (38) 

and we have the following properties: 

H(txM)= f dtF(ti,t)G(t,t2)^H(a;i,a;2) = J]^(^i''^)^("^'^2) (39) 



• \ 2 



H{tiM) = F{tiM)G{tut2)^n{uu<^2)= Y{uj[,uj'^)G{uji-uj[,uj2-uj'^) (40) 



The transformed free scalar propagator is: 



Do(^l,^2) = -5uj^+u,^,0 X2 _ 2 (41) 



All these relations enable us to use the usual "ITF" Feynman rules. The important point 
to note here is the fact that these Feynman rules are not at all dependant of the choice of the 
path going from ti to ti — i/3, so that the name "Imaginary Time Formalism" gives a very 
restrictive view of this formalism. 
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4.4 Conclusion 



The conclusion of this part is that every graph evaluated on the complete path C is invariant 
by time translation and does not depend explicitly on the initial time tj. In the proof, this 
property has been shown to be closely related with thermal equilibrium through the KMS 
relations. As a corollary of this result, we can look with suspicion the arguments based on the 
limits ti — oo, tp +00 because nothing special occurs when we look at these limits. 



V. The vertical part contribution for two point functions 



5.1 Spectral representation of two points functions 

In order to extend the result of part III to general two point functions, we use their spectral 
representation. Let ^(ii, ^2; k) be a two point function; we can write it as [3 ,10] : 

+00 

HtiM-.k)= j dE f{E,k)Goiti,h;E) (42) 



where f{E, k) is its spectral density and Gq{-) the free propagator. In particular, this is true 
for a self-energy function: 

+00 

^{h,t2;k)= J dEa{E,k)Go{ti,t2;E) (43) 


Then, if Ejj is the matrix self-energy, we have: 

+ CXD 



Sii(A;o,fe)= J dE a{E,k)Du{ko,E) ^22{ko,k)= J dE a{E, k)D22{ko, E) 



+00 +00 

Si2(fco,fe) = - J dEa{E,k)Di2{ko,E) ^2i{ko,k) = - J dE a{E,k)D2i{ko,E) 



where the minus signs in the 12 and 21 components are due to the opposite sign of the type 2 
vertices. 
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5.2 The insertion of a self-energy term 
5.2.1 Calculation in time coordinates 

We now consider the insertion of a self-energy type function between two free propagators^"' 
(figure 4): 

F{ti,t2,k) = j dtzdUGQ{ti,U-ujk){-i^{tzM\k))GQ{UM\^k) (45) 



Thanks to the spectral representation of S(-), this calculation is simply the convolution of 
three free propagators. In the following, we need only the contribution of the vertical part, 
which is: 

+ 00 

Its Fourier transformation is independant of the ij component considered because Fy(-) does 
not contain any 9c{-) function. We have: 

^.„,.^.(,„, = |!!£W(1±!!£M)2.*(*S - .1) J iE .(E, k)^ (47) 



5.2.2 Calculation in RTF 

Now, we do the same calculation with the usual Feynman rules of RTF in order to see the 
difference between the nsiu^k) and the ns(|/eo|) prescriptions. We have to evaluate: 

Fu,[ko],ijiko, k) = ^ D^iko],iaiko,uJk)i-iT,abiko,k))D^iko],bjiko,i^k) (48) 

a, 6=1,2 

If we explicit a bit more this quantity, we see that the ill defined products of distributions 
cancels mutually thanks to the following identities: 

fSii + E22 + S12 + S21 = 
ie{-ko) + nBi\kom2iiko,k) = {e{ko) + nB{\ko\))^i2{ko,k) ^ ^ 

(The first one is always true, and the second one, being a form of the KMS relation, is only 
true for equilibrium). 

Moreover, this cancellation makes derivatives of 5{-) functions appear. For exemple, the 11 
component is: 

iV 

Fco[ko]Mko,k) = —2 -2N2^1l(^0'^) 

^ ^ (50) 
- TT [(1 + 2nB)En + {e{-ko) + nB)^2i + (^(fco) + nB)^i2] -^K^l - ool) 

OKq 



This operation is the basis of the construction of the complete propagator. 
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Therefore, the difference Fk^^^ij{ko, k) — F^^ij{ko, k) will not be zero (see appendix 3). By 
an explicit calculation, we verify that this difference is independant of the component ij one 
is considering: 



Fko,ijiko, k) - F^^ij{ko, k) = 27rS{ko - w^; 



5.3 Epilogue 



-i En + 



(51) 



Then, thanks to the spectral representation of the self-energy matrix and by using the 
explicit form of the free matrix propagator (see appendix 2), we can verify the identity: 



+ 00 



Sn + ^ = ^ = J dE a{E, k)-^ (52) 



So that we have: 

Fy^ij{ko, k) = Fko,ij{ko, k) - F^,ij(fco, k) (53) 

So, we have proven that the contribution of the vertical part corresponds exactly to the dif- 
ference between the results given by the nB{\ko\) and the nB{uJk) prescriptions. 

Therefore, to take this vertical part into account, we have to apply the n^dfeol) prescription 
in the usual Real Time Feynman rules'^°\ at least for two point functions. 



VI. General remarks on n point graphs 
6.1 Preliminaries 

Before going to the case of higher n point functions, let's consider the following integral'^^': 

jr(E) = J dtF{t)e-'^^ (54) 
c 

where E G R and where F{-) is the characteristic function of a subset ^ of C {F{t) = if 
t ^ A, F{t) = 1 if t e ^), with Af] Cv = or C„. We will explain later the reason of this 
condition. Our purpose is here to determine the contribution of the vertical part in If 

10) Here, we admit that the contribution of the two horizontal parts is obtained by the usual 
RTF Feynman rules with 77.^(0;^.) prescription. This point is rather obvious because the prop- 
agators in time coordinates already contains nB{u}k) distribution functions. 

This is the kind of calculation we have to perform in time coordinates, if we except the 
three momentum integrations in loops. 
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Ad Cy = 0, this contribution is trivially zero, so that we only consider the other case. We 
have now to distinguish between E ^ and E = 0: 

(i) if E 7^ 0, we have: 

JF(E) = ^ e-^^*Me-^^ - 1) + J dtF(t)e-'^' (55) 

But we know that the result of the calculation of the graph, which is a sum of such terms, 
is independent of tj. So, there must exist other terms which cancel the contribution of the 
vertical part in (55), because it depends on tj. The conclusion is that when E 7^ 0, (54) will 
not give a vertical part contribution to the graph. 

(ii) if E = 0, we have: 

jr(E) = -i/3+ J dtF{t)e-'^* (56) 

Now, the preceding argument cannot be applied because the contribution of the vertical part 
does not depend on tj. So, in general, there is a contribution of the vertical part in (54) to 
the graph when S = 0. 

6.2 Application to n point functions 

The calculation of the vertical part contribution to the time dependence of a n point function 
can be reduced to a sum of integrals like (54), where F(-) is a product of 9c{-) functions and 
E a linear combination of the energies cui = •>/ m? + k? carried by the various lines of the 
graph. The hypothesis made above on ^ fl Cy is related to the fact that we have defined the 
contribution of the vertical part to be the one obtained when the integration time is on C^, 
and the other times on Ci © C2 (see note 5). 

So, we have shown that the result of a vertex integration yields a vertical part contribution 
if and only if its E is zero. Moreover, this contribution is always finite {—if3 midtiplied by 
a regular function of the other times). Then, if we take into account the fact that we have 
to perform integrations over the internal independant three momenta, this contribution will 
vanish at the end when E contains an internal momentum because E = defines a null mesure 
subset of the set in which this momentum can vary. 

The only remaining possibility to have contributions of the vertical part is then in terms 
for which E contains only the external three momenta, or if the diagram contains a self-energy 
subdiagram (this situation has been studied in the preceding section). 

6.3 Example 

In order to illustrate this, we consider the example of figure 5. The external propagators, 
carrying energy uji, are denoted for example: 

Go(ti,«i;wi) = ^ J2 [^c(ei(«i-ti)) + nB(a;i)]e-™(-^-*^) (57) 

ei=±l 
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and the internal ones are denoted for example: 



So, we have to perform an integral like: 

JflidUiC-'^'''') F{ti,t2,t3;Ui,U2,U3,U4,U5) (59) 



c 



i=l 



with: 

{El = eiWi - ?7ifii + 776^^6 , S2 = €2U!2 - ri2^2 + Vl^l ; ^3 = ^3i^3 " Va^S + 112^2 ^^^^ 

At the first integration (over mi), the S introduced above is Si. At the second integration, it 
can be S2 or Ei + E2 because the first integration can generate a term in exp(— zEitt2) . At 
the third integration, it can be E3, Ei + E3, E2 + E3 or Ei + E2 + E3, and so on. Among 
all the possible E encountered in the five integrations, only one depends only on the external 
energies tUi: 

3 

El + E2 + E3 + E4 + E5 = J] CiiOi (61) 

i=l 

which is encountered in the last integration. Moreover, because of the conservation of three 
momentum, we have fi3 = fie so that: 

E4 + E5 = (7/3 - ve)^3 (62) 

and can be zero independantly of fl^ when 773 = 776: this is precisely what occurs when a 
subdiagram is of self-energy type. 



6.4 Conclusion 

The conclusion of this section is that it is possible to have a contribution of the vertical 
part of the path in > 2 point functions for certain configurations of the external momenta 
(namely, when a combination Yli=i ^i^i of the external energies is zero) or when the graph 
contains a self-energy subdiagram. 
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VII. Concluding remarks 



The main result of this paper is to exphcit the role of the vertical part of the real time path, 
giving a rather physical interpretation of its role. We proved that the vertical part is necessary 
in order to have the cancellation of the terms which are not invariant by time translation, 
and that the remaining terms do not depend on the initial time choosen for the path. We 
must emphasize that such a property was to be expected on physical considerations for an 
equilibrium system. As a consequence, the usual proofs of the Fourier Feynman rules based 
on the limits tj — > — oo, tp — > +00 are thought to be false. 

The other result is that the vertical part of the path generally contributes to the Green 
functions, as already seen in [10]. Moreover, we have shown in the case of two points functions 
that this contribution is simply taken into account by the nB{\ko\) prescription associated to 
the usual Real Time Feynman rules. This proof is consistent with the fact that the Green 
functions does not depend on the contour parameters, because it works even if we keep ti 
finite. 

It remains to study more precisely the case of the functions with n > 2 points, and to see 
if the contribution of the vertical part which can occur can be dealt simply by the nB{\ko\) 
prescription. Another problem is to make more precise the relation between pathologies and 
the fact that the theory is or is not an equilibrium one, as outlined in [17,18] . 
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Note added in proof: 

After this work had been completed, we became aware of a preprint by T.S. Evans and 
A.C. Pearson [19], which contains a rather different approach to the problem of the vertical 
parts of the path in RTF. 

In [19], an e— regularized propagator is used: 

Gt,{h,t2;L0k) = Go(ti,t2;a;fc)e-^l*^-*il (63) 

which vanishes when the time difference goes to infinity, in association with the limit tj — > —00. 
In order to clarify the contradiction between their result and ours concerning the contribution 
of the vertical part, we repeated the simple calculation presented in section III of the present 
paper with their regularized propagator, and we were led to the conclusion that the limits 
e — > 0"*" and t/ ^ — 00 do not commute because of terms like exp(et7-). More precisely: 

(i) If we do e ^ 0+ first (this is equivalent to the use of non-regularized propagators) , the 
result is invariant by time translation and does not depend anymore on tj, so that the limit 
t/ ^ — 00 is useless. Moreover, there remains a contribution of the vertical part of the path in 
the result. 

(ii) If we do t/ ^ —00 first, like in [19], the result is also invariant by time translation (but 
this time for different reasons, because the free regularized propagator does not verify KMS) 
but the contribution of the vertical part is lost. 
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Because of the physical interpretation of thermal equilibrium, we believe that in a consistent 
answer to this problem the limit tj ^ —oo has no role to play, so that in time coordinates, 
the use of non-regularized propagators seems preferable. 



Appendix 1 : The generating functional 

Our purpose is to transform expressions like: 

Z^'ij] =uj ^Tccxpz j dt^x Ci{Mx))+jix)(t)j{x) 

= exp i / d'^x Cj 



7 

d 



idj{x) 



7 



To achieve that, we will use twice the Campbell-Haussdorf formula in order to transform the 
path ordering in normal ordering [14]. By doing this, we make a c- number appear, and the 
main point is that a c-number is equal to its average c = uii{c), where a;i(-) is another arbitrary 
average operator. The general result is then (in the case of scalar fields): 

Z-[j] =expij d'x Cj {Jj^) X Z^IJ] iA2) 



with 



Z^j] = C[j] X exp-i J d^xd^y jix)j{y)LUi{T,Mx)My)) (^3) 

7 

X exp^ J d'^xdS j{x)j{y)uji ( (^+(x)07(y) + 0+(y)07(x) 



(AA) 



(j>+{x)(l)+{y) + 4>j{x)4>j{y) ) 



where (pf and (jyy are respectively the creation and the annihilation part of the free field^^'. 
The first factor in A2 will be expanded in powers of the coupling constant in order to perform 
the perturbative expansion. The factor C[j] contains the information relative to the initial 
correlations of the system: we have C[j] = 1 when the Wick's theorem holds. The equilibrium 
field theory is a bit particular because we can replace ujh{-) by ujho{-) in C[j] provided we 



This expression for ZqIJ] is consistent with the Dyson- Schwinger equation for real scalar 
fields [3]: 



.2jZ-[j] 



- 17 - 



add the vertical part Cy to the path Ci © C2. Then, by choosing coi = uJhq'^^\ we can make 
the correlation terms disappear'^' , so that the Feynman rules will be the usual ones, the only 
difference being the integration path. 

We can add that in the general case, there is no simple choice of the path and of the 
auxiliary average allowing C\j] = 1. The consequence of this impossibility is that one will 
have very intricated Feynman rules with not only propagators, but also correlation functions 
[20,21,22] . 

Appendix 2 : The free RTF propagators 

Here we recall the free matricial propagators we use in RTF with a = 0: 

iP 

-Da,[feo], 11(^0, t^fe) = Ta 2 +'^i'^ + '^^B{ojk[\kQ\]))5{kl- u)l) 

—iV 
'^0 ^k 

Duj[ko],i2{kQ,ujk) = 27r{e{-ko)+nB{uJk[\ko\]))5{kQ-ujl) 
,-D^[feo], 21(^0, t^fc) = 2'K{e{ko) + nBiuJk[\ko\]))5{kl - u^l) 

where the argument of the distribution functions is cuk or \ko\, according to the prescription 
we choose. 

Appendix 3 : Distributions and pathologies 

The distributions we have to deal with in this paper are S{x) and V/x. The squares 5^{x) 
and (V/x)'^ are ill defined, even in terms of distributions, but it is possible to give a sense to 
a linear combination of them. To do so, the proper method is to regularize the distributions 
we are interested in. A simple regularization is given by [3]: 

lim 

where the derivative is defined in the sense of distributions: 

V/, D'if) = -Dif) 
By using (A5) for n = and n = 1, we obtain the following relations: 

^=(^y-n'5\x) m 
\x J 

The fact that we can make C[j] = 1 with loi = cuHa is very important to have a consistent 
free theory. Indeed, ujHaC^c4'i{x)4'i{y)) is a free quantity both for the dynamics of the fields 
and for the density operator with which the average is performed. 

The denomination "correlation terms" is a bit ambiguous in this context because they 
depend on the the actual choice of the integration path and of the auxiliary average 
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{A5) 



X + ie 



X 



V 

n+l 



^ .7r(__y 



(A6) 



iA7) 



26{x)- = -5'{x) (v49) 

X 

These relations, combined to the definition of the distributional derivative, imply that the 
prescription for the argument in the statistical factors is important each time we have to 
multiply free propagators carrying the same momentum. Indeed, if we look at a product like 
a{x)5{x){P / x), the result will be a different distribution if we replace 2; by in a(-)^'*'. In order 
to see that*^', consider the distribution {a{x) — a{Q))5'{x) and apply it to a test function /(x); 
the result is not zero: 



+00 +00 

dx (a(x) - a{Q))5\x)f{x) = - j dx 5{x)[{a{x) - a(0))/(x)]' 
00 —00 

= -a'(0)/(0)7^0 



Of course, we have a{x)5{x) = a{0)5{x). In fact, the problem comes from the fact that we 
cannot define unambiguously an associative multiplication of distributions. 
'•^^ Another way to see the problem is to perform the inverse Fourier transformation of a 
function where pathologies does arise. Then, one will encounter double poles, the residue 
of which contains derivatives: the various prescriptions for the arguments of the distribution 
functions will lead us to different results. 
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FIGURE CAPTIONS 



Figure 1 The Real Time Path 

Figure 2 The Real Time Path with a = 

Figure 3 Repeated insertions of mass terms 

Figure 4 Insertion of a self-energy function 

Figure 5 An example of 3 point diagram 
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